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We investigate the bipolaronic crossover and the pairing transition for a two-orbital model with Jahn-Teller 
coupling to a two-fold degenerate phonon mode. The evolution from weak to strong coupling is reminiscent 
of the behavior of the single-band Holstein model: The polaron crossover in which the electrons and phonons 
become strongly entangled occurs for a weaker coupling than the binding of bipolarons, which gives rise to 
a metal-insulator transition. Interestingly, a single bipolaronic transition takes place also when the two bands 
have significantly different bandwidths, as opposed to the case of repulsive Hubbard-like interactions for which 
an orbital-selective Mott transition has been reported. This behavior is related to the inter-orbital nature of the 
Jahn-Teller coupling. 
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I. INTRODUCTION 

The interplay between spin, charge and orbital degrees of 
freedom plays an essential role in the description of the com- 
plex phase diagrams characterizing strongly correlated mate- 
rials. One of the important mechanisms describing the prop- 
erties of these systems is the electron-phonon (e-ph) inter- 
actions. Most of the studies of e-ph interactions are limited 
to the Holstein model, in which the electronic state forms a 
non-degenerate band and the phonons couple with the on-site 
electronic charge 1 . However, many strongly correlated mate- 
rials involve nearly degenerate d or / orbitals. In this work 
we move a step towards a more realistic description of the 
e-ph interaction in these compounds by considering a degen- 
erate electronic manifold coupled with the proper degenerate 
phononic models. For these narrow-band systems the phonon 
dynamics can play a role, therefore the Born-Oppenheimer 
approximation is hardly justified and the decoupling of the 
electronic orbital states from the vibrational modes is not al- 
lowed. 

A typical family of compounds involving such a spin- 
charge-orbital complex are the colossal magnetoresistance 
(CMR) manganites 2 . While the double-exchange mechanism 
is the basis of the CMR phenomenon, relating the magnetic 
behavior to conduction properties, its interplay with the cou- 
pling between the degenerate e g electrons and the Jahn-Teller 
(JT) distortions of the MnOg octahedra is crucial to describe 
the properties of these materials 3 leading to a variety of exper- 
imentally observed charge and/or orbital orders 4 . 

The Jahn-Teller model is also relevant to superconductiv- 
ity in alkali-doped A^Cgo molecular solid, where Ceo is the 
fullerene molecule and A stands for alkali cations K, Rb, 
or Cs. In the alkali-doped A x C6o, the threefold degenerate 
t\ u molecular level is partly occupied and couples strongly 
to eight H g intra-molecular Jahn-Teller phonons 5 . In these 
systems the JT character of relevant phonons has important 



consequences in the presence of strong Coulomb repulsion 6,7 . 
Indeed the JT-driven interaction between the electrons only 
involves spin and orbital degrees of freedom which are still 
free to fluctuate when the charge excitations are frozen by the 
Coulomb repulsion. This gives rise to a correlation-driven en- 
hancement of phonon-driven superconductivity 6 . 

In both cases the interplay between electron-electron and 
e-ph interactions which may lead to a rich variety of physical 
phenomena. In this paper we focus on the pure e-ph inter- 
action term in order to highlight the specific properties intro- 
duced by orbital degeneracy before considering the explicit 
role of electronic interactions, in the same spirit of previous 
investigations of the Holstein model 8 ' 910 . For the same rea- 
son, i.e., capturing the basic aspects of the electron-phonon 
interaction we do not allow for superconductivity or charge- 
density-wave ordering, limiting ourselves to the normal state. 
In the half-filled Holstein model it has indeed been shown that 
two related but distinct processes occur by increasing the e-ph 
coupling. The first effect is the well-known polaron forma- 
tion, i.e., the progressive entanglement between the electronic 
motion and the lattice degrees of freedom. Polaron forma- 
tion occurs indeed as a continuous crossover and it can be 
pinpointed by analyzing the phonon displacement distribution 
function. The attractive interactions also induce an attraction 
between the fermions. Even if superconductivity is explic- 
itly inhibited, this gives rise to a binding of fermionic carriers, 
eventually leading to a pairing transition 12 which forpolaronic 
carriers becomes a bipolaronic metal-insulator transition 8,9 . 

In this paper we discuss how this physics is modified for a 
two-orbital e x E Jahn-Teller model. For this model we can 
also consider some perturbations, like the possibility of differ- 
ent bare widths for the two electronic bands, in analogy with 
the much discussed possibility of an orbital-selective Mott 
transition in the case of repulsive Hubbard-like interactions 13 . 

The paper is organized as follows. After briefly introduc- 
ing the model and the method of solution in the next section 
(Sec. II), we begin to investigate the JT model at half-filling 
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in Sec. Ill and discuss how the e-ph coupling affects the elec- 
tronic properties of the model. In Sec. IV, we then address the 
question of whether orbital-selective bipolaron transition oc- 
cur in the system with different band-width. The concluding 
remarks is given in Sec. V. 

II. MODEL AND METHOD 

The e x E Jahn-Teller model consists of two degenerate 
electron orbitals and two degenerate phonon modes, 

H = H t + H p h + Hjt 

(ij)ja- 

H ph = Q y^qtqj + b\bj) 

i 

Hj T = gy^\{nn a - n, l2a ){a\ + a,) + 

icr 

+ {c\ la c aa + c\ 2a c lXa ){b\ + &;)], 

(1) 

where a\ (bj) and a>i (bi), respectively, create and annihilate 
the dispersionless phonons of type a (b) on site i. These are 
the Jahn-Teller modes with the same symmetry and the same 
phonon frequency i7o- g is the related electron-phonon cou- 
pling strength. While phonons of type a couple to the density 
unbalance between the two orbitals, the phonons of type b are 
coupled to the hybridization between electronic orbitals that 
are orthogonal in the absence of these phonons 14 . c| 7(T and 
Ci 1<y are the creation and annihilation operators for electrons 
at site i in orbital 7 (7 = 1,2) with spin tr, n ila = c\ ia a lcr 
and (ij) denotes nearest neighbors. t 1 > are the nearest- 
neighbor hopping integrals. 

We solve the model by means of dynamical mean field 
theory (DMFT) 15 . The method maps the lattice model onto 
a quantum impurity model in which an interacting site is 
embedded into a non-interacting bath, whose spectral func- 
tion has to be determined self-consistently in such a way 
that the impurity Green's function of the quantum impurity 
model coincides with the local Green's function of the lat- 
tice model under consideration in order to enforce a quan- 
tum dynamical mean-field theory. This requirement leads 
to a self-consistency condition, which contains the informa- 
tion about the original lattice through the non-interacting den- 
sity of states. A particularly popular and useful choice is 
a semi-circular density of states of half-bandwidth D, i.e., 
N(u>) = -^jyi^/D 2 — oj 2 , which corresponds to an infinite co- 
ordination Bethe lattice. For this system the self-consistency 
equation takes the following simple form: 

r>2 V 2 

— G^„)=^-^, (2) 

where G(iui n ) is the local Green's function of the system, 
and Vfe are the energies and the hybridization parameters of 



the impurity model. To solve the impurity model at zero tem- 
perature (T = 0) we use the exact diagonalization method 16 , 
which works equally well at any value of the parameter and 
only involves a discretization of the bath function, which is 
described in terms of a finite and small set of levels n s (con- 
ventionally n s includes also the impurity site and the number 
of bath sites is rib = n s — l) in order to limit the Hilbert space 
to a solvable size. 

In principle, an infinite number of vibrons can be excited 
at each lattice site. In DMFT, however, the phonon degrees 
of freedom are limited to the impurity site. Since the diag- 
onalization can only be carried out for a finite dimension, a 
truncation of the local phonon subspace is required. Here we 
use the basis state 

W) P h ; v v = I y j n it u < N ph (3) 

leading to (N ph + N mode )\/{N ph \N mode \) allowed phonon 
configurations. Here m :V is the number of i-type phonon in 
the basis state \v) p h and N moc i e = 2. Typical values we con- 
sidered for the bath levels are n s = 5 — 6 (which corresponds 
to n s = 10 — 12 for a single-band model) and N p h ~ 30 — 40. 
We tested that these numbers provide essentially converged 
results. 

DMFT has already been applied to the study of strongly 
correlated electron-phonon systems and has emerged as one of 
the most reliable tools for the analysis of these systems 8 ' 9 ' 10 ' 11 . 
Studies on the normal phase of the Holstein model show that 
for e-ph couplings, the ground state is metallic with Fermi liq- 
uid characteristic, while increasing the coupling a first-order 
metal-insulator transition takes place. This phase displays a 
gap in the one-electron spectrum, but, in contrast to metal- 
insulator Mott transition, no hysteresis and no preformed gap 
is observed for the transition from the metallic to the bipola- 
ronic insulating phase, at least for the relatively small values 
of phonon frequency that have been used. The pairing tran- 
sition is accompanied by, but it does not coincide with, the 
progressive entanglement between electronic and vibronic de- 
grees of freedom, or in other words the polaron crossover. For 
relatively small phonon frequencies, polarons are formed be- 
fore the pairing transition occurs. Since the latter takes tran- 
sition binds two polarons, it can be defined as a bipolaronic 
transition 8 . One of the aims of this paper is to understand 
whether and to what extent the symmetry of the e-ph coupling 
can influence the electronic properties and the dynamic quan- 
tities of a given system. 

To analyze the metal-insulator transition we study the 
quasi-particle weights Zi = l/[l — dTZeSi(u>) / du)\ u =o\ > 
being the self-energy for band i. Due to the momentum in- 
dependence of the self-energy in DMFT, z is proportional to 
the inverse effective mass (z oc m/m*), so that a vanishing z 
implies a divergent effective mass and a metal-insulator tran- 
sition. 

We will also use the double occupancy dj = (n^riii), 
and the phonon displacement probability distribution function 
(PDF), defined by: 

P(x) = (ip \x}(x\i> ) , (4) 
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as a marker of the polaron crossover 8 ' 9 . In Eq.(4), |^q) is 
the ground state vector, and |a;)(x| is the projection oper- 
ator on the subspace where the phonon displacement value 
x has a given value x. This quantity, therefore, is a mea- 
sure of the distribution of the local distortions. In the 
exact diagonalization approach we can evaluate P(x) = 
Sn,m 0n 0*0 0m (x) (-00 1 n) (m\ ipo), where \n) and |m) are the 
eigenstates of the harmonic oscillator and \(f> n (x)) are the cor- 
responding eigenfunctions. 



III. RESULTS 

We first focus on the isotropic systems consisting of identi- 
cal bands, D\ = D2 = D. Furthermore, we confine ourselves 
to the adiabatic regime with a small ratio of the phonon fre- 
quency to the semi-bandwidth (Oq = 0.2Di) in which the 
dynamics of phonons is slow as compared to the typical ki- 
netic energy of electrons. The results obtained for this case at 
half-filling are shown in Figs. (1) and (2). In Fig. (1), panel 
(a) we report the quasi-particle weight z; = z as a function 
of electron-phonon coupling constant. In a system with e-ph 
interaction, by increasing the coupling constant, the system 
enters a strong-coupling polaronic regime in which the pres- 
ence of an electron is associated with a finite lattice distortion. 
In the polaronic regime, the system is characterized by low- 
energy excitations which, although strongly renormalized, are 
still coherent. Also, the same e-ph coupling can cause a pair 
of polarons with opposite spins to attract each other forming 
a bound state in real space, called bipolaron. In DMFT those 
pairs are unable to move, and the bipolaron formation causes 
a metal-insulator transition associated with the divergence of 
the effective mass and the disappearance of a coherent peak 
at the Fermi level at a critical coupling g c . This process is 
evident in our DMFT results from the enhancement of the 
electron effective masses, and the corresponding reduction of 
the quasi-particle weights (zi = Z2 = z) as the e-ph cou- 
pling constant g increases, which eventually reaches the bipo- 
laronic metal-insulator transition (MIT) shown in panel (a) of 
Fig. (1). Panel(b) of Fig. (1) shows the double occupancies 
{di = d 2 ) as a function of g. Around the critical coupling, 
g c , the double occupancies increase suddenly, showing a first 
order transition to bipolaronic phase. 
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FIG. 1: Panel (a): The quasi-particle weights, zi = 22, at half-filling 
as a function of electron-phonon interaction, g. Panel (b): double 
occupancy, di — di, as a function of g. The inset of the panel (b) 
shows the orbital correlation, (runs), between the two orbitals. 



The inset of the Fig. (1), panel (b), shows orbital corre- 
lations between the two orbitals, (711712). Upon increasing 
the e-ph coupling, the orbital correlation ((ni}(n2)-(nin2) = 
1 — (7)1712)) increases until bipolaron transition takes place. 

Evolution of quasi-particle weights and the double occu- 
pancies do not clearly show the polaron crossover. There- 
fore, we have computed the phonon displacement probability 
distribution function to clarify this point. Fig. (2) show the 
evolution of the phonon displacement probability distribution 
function, P(x), as a function of g. In the case of isotropic 
systems, the PDFs for a- and b-modes are the same. Upon 
increasing the e-ph coupling, a smooth crossover occurs be- 
tween a unimodal distribution and a bimodal distribution in 
complete analogy to the single-band Holstein model 8,9 . 




g-0.00 
g.0.10 

g.0.14 
g.0.18 
g.0.20 
g-0.232 
g-0.238 
g.0.239 



6.0 -4.0 -2.0 0.0 2.0 4.0 6.0 



FIG. 2: Phonon displacement distribution function P(x). The vari- 
ous lines correspond to different values of g. 



Before entering the bipolaronic phase characterized by z = 
0, a very narrow region exists in which polarons are formed, 
nonetheless, the value of quasi-particle weights, although be- 
ing very small, are different from zero (see two dot-dashed 
curve of Fig. (2)). The existence of this region shows that, in 
the process of the metal-insulator bipolaronic transition, bipo- 
larons are formed for those values of g which are somewhat 
larger than the polaron crossover. While the qualitative result 
is similar to what found for the Holstein model, the region in 
which polarons exist without forming bipolarons turns out to 
be narrower in the Jahn-Teller model 8,9 . 

We have so far discussed the bipolaronic transition for the 
two-orbital electronic system with the same bandwidth. How- 
ever, as a result of the geometric complexity of many tran- 
sition metal oxides, the degeneracy of the valence bands is 
frequently lifted, giving rise to a coexisting partially filled 
narrow- and wide-bands. Actually, the two-band Hubbard 
model with orbitals of different widths has recently received 
considerable attention and it has been shown that in certain 
range of repulsive Coulomb interactions and Hund's cou- 
pling the so-called orbital-selective phase occurs in which 
the narrow-band is insulating while the wide-band is still 
metallic 13 . 

In the case of a phonon interaction involving the total 
charge on the two orbitals, as in a generalized Holstein model, 
the situation would be most likely similar to the case of the re- 
pulsive Hubbard model. Indeed a Holstein coupling induces 
a charge-charge attraction which, in the antiadiabatic limit 
Slo — ^ oo reduces to an attractive Hubbard model, which in 
turn is equivalent to the repulsive one at half-filling 12 . For the 
JT interaction the question is instead not obvious because the 
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phonon modes are not coupled to the total charge, and they are 
associated with the orbital degrees of freedom. In particular, 
the 6-mode is directly associated with a hybridization between 
the two orbitals, an effect which is clearly against an orbital 
selective behavior. If we integrate out the phonons in our ef- 
fective action in which the interaction part of (1) is present 
only on the impurity, we obtain a retarded attraction 



V eff {u) = g 2 



+ 



E 



2>M 



(5) 



where T>{lo) = — 2£lo/(^o — uj2 ) 15 ^ e ^ are phonon propaga- 
tor for each mode. In the antiadiabatic limit fio — * oo we can 
drop the frequency dependence in the propagator obtaining an 
unretarded attraction which involves the orbital momentum of 
the electrons in the two orbitals 23,24 and coincides with an in- 
verted Hund's rule term 6 . 




FIG. 3: Quasi-particle weights, zi, 22, (left panels) and double 
occupancies, di, eta, (right panels) at half-filling as a function of 
electron-phonon interaction for different bandwidths. Top panels: 
D2/-D1 = 0.5, bottom panels D2/D1 = 0.2. The insets of right 
panels show the orbital correlations, (rixriz) as a function of g. 



with respect to their bare values. However, when the sys- 
tem approaches the polaron crossover, the difference between 
quasi-particle weights reduces, displaying a fast damping de- 
pendence on g. Therefore we have a single bipolaronic transi- 
tion, at least within our numerical accuracy in a wide range of 
D2/D1 and the orbital-selective effects are extremely limited 
if compared with charge interactions. This is interesting also 
because reducing the value of D2, the second band is not in 
the same adiabatic regime of the broad one, which could imply 
some different behavior. Nonetheless, the inter-orbital nature 
of the e-ph coupling overcomes also this effect to leading to a 
unique transition. The spectral functions of the two bands (not 
shown) confirm the uniqueness of the transition. At the bipo- 
laronic transition, where the two z's vanish simultaneously, 
the peak at the Fermi level disappears in both bands. 

On the other hand, the growth in the number of doubly oc- 
cupies sites (right panels in Fig. 3) depends on the band and 
on the value of the bandwidth ratio. The narrow band rapidly 
develops a larger number of doubly occupied sites, and the 
effect is more pronounced the narrower the band is. The crit- 
ical value of g for the metal insulator transition decreases as 
the ratio D2/D1 is reduced. This behavior is somewhat nat- 
ural since we are keeping D\ fixed and reducing £>2- There- 
fore the overall kinetic energy of the system is reduced, mak- 
ing polaron and bipolaron formation easier (both the polaron 
crossover and the bipolaron transition are expected to occur 
because the e-ph interaction overcomes the kinetic energy and 
its delocalizing effect, at least in the adiabatic regime in which 
the energetic convenience also implies a sizeable entangle- 
ment between electrons and phonons 17 ). 

The phonon displacement distribution functions are shown 
at Fig. (4) for D2/D1 = 0.5. Breaking the symmetry be- 
tween the two orbitals determines a different behavior of the 
two phonon modes. With increasing e-ph coupling, the PDF 
of the a mode becomes bimodal before the bipolaronic tran- 
sition occurs, while for the same values of the coupling the 
PDF of the 6-mode maintains a unimodal shape with a single 
feature which broadens without showing bimodality all the 
way to the bipolaron transition. The two peaks of the bimodal 
distribution for the a mode turn out to have slightly differ- 
ent height. The width of the region of couplings between the 
polaron crossover for the a mode and the metal insulator tran- 
sition is larger than for the case of equal bandwidths. 



In order to address the existence of the selective bipolaronic 
phase, we considered two different ratios between the bare 
bandwidths: D 2 /Di = 0.5 and 0.2. Fig. 3 shows the re- 
sults for quasi-particle weights and double occupancies of the 
two bands in the two cases (The top panel is D2/D1 = 0.5 
and the bottom one D2/-D1 = 0.2). The results for the z's 
are apparently similar to those shown in Fig. 1 for identi- 
cal bandwidths. The two renormalization factors are indeed 
only slightly different, in sharp contrast with the results for 
the Hubbard model 13 . Moreover, the difference in renormal- 
ization goes in the opposite direction than for repulsive mod- 
els with pure charge interactions. Here the wider band is 
(slightly) more renormalized by the interactions, as if the ef- 
fect of e-ph coupling is to make the two bandwidths closer 
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FIG. 4: Phonon displacement distribution function P(x). The left 
panel is for a phonon mode which shows bimodality. The right panel 
is for b phonon mode which does not show bimodality. 
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IV. CONCLUDING REMARKS 

We considered the normal phase of a two-orbital model 
coupled with a two-fold degenerate phonon manifold via an 
e x E Jahn-Teller coupling. Solving the model at half-filling 
and zero temperature by means of dynamical mean-field the- 
ory with exact diagonalization as the "impurity solver", we 
investigated the polaron crossover and the bipolaron transi- 
tion which occur as a function of the coupling constant. In the 
standard case of identical bands, our results are qualitatively 
similar to those of a single-band Holstein model. In the adi- 
abatic regime that we investigate, a polaron crossover takes 
place for weaker coupling with respect to a metal-insulator 
transition associated to the formation of bipolarons 89 . The 
intermediate region between the two effect is found to be 
smaller than in the single-band case. 

Interestingly, the situation is essentially unaltered if the two 
bands have different bandwidths. The system undergoes a sin- 
gle bipolaronic transition, which is preceded by a polaronic 
crossover for one of the two phononic modes. 

Nonetheless, the effect of electron-electron interaction can 
be significantly different in the present two-orbital case. In 



the single-band case, a charge-charge Hubbard repulsion di- 
rectly competes with the Holstein term which couples the 
phonons to the local charge. This gives rise to a phase diagram 
characterized by the competition between phonon-driven and 
correlation-driven localization effects 18 ' 19 ' 20 ' 21 ' 22 . 

In the present case, the phonons couple to electronic de- 
grees of freedom other than the total onsite charge. This 
makes the competition between the two terms much more sub- 
tle. In a purely electronic model which coincides with (1) 
if the antiadiabatic limit is taken for phonons it has been in- 
deed shown that phonon-mediated superconductivity can be 
enhanced by correlations 23 and several anomalous properties 
can arise 24 . The investigation of the combined effects of cor- 
relation and JT coupling in the present model is therefore ex- 
pected to lead to a rich and interesting physics. 
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